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\HsiRAt i: The potent dipeptidyl peptidase IV (DP IV) inhibitor 1 1 -( 2-py rrolidiny Icarbonyl )-2-py rroli- 
din\l|boronic acid (t.-Pro-Di -boroPro) (I-lcntkc. G R , Munoz. E , Hubcr, B. T., Plaut, A G , Kettner, 
( A., & Bachovchin. W. W. (1991) Proc. Sail. Acad. Sa U.S.A. SH. 1556- 1559] was fractionated into 
its component l-l and l-d diastereomers by C 1 8 HPLC, and the binding of the purified diastereomers to 
DP IV was analyzed. Inhibition kinetics confirms that the l-i diastcreomer is a potent inhibitor of DP IV, 
having a A.', of 16 pM The L-D isomer binds at least 1000-fold more weakly than the L-L, if it binds at all, 
as the - 200-fold weaker inhibition observed for the purified l-d isomer is shown here to be due entirely 
to the presence cf a small amount (0.59 p /r ) of the l-l diastcreomer contaminating the l-d preparation. The 
instability of Pro-boroPrc. together with its very high affinity for DP IV and the time dependence of the 
inhibition, makes a rigorous kinetic analysis of its binding to DP IV difficult. Here we have developed a 
method which takes advantage of the slow rate at which the inhibitor dissociates from the enzyme. The 
method involves preincubating the enzyme and the inhibitor without substrate and then assaying the free 
enzyme hy the addition of substrate and following its hydrolysis for a period of time which is short relative 
to the dissocia: on rate cf the inhibitor. Data from experiments in which the preincubation time was 
sufficient for enzyme and inhibitor to reach equilibrium were analyzed by fitting to an appropriate form 
of the quadratic equation and yielded a K, value of 16 pM. Data from experiments in which the incubation 
lime was insufficient to establish equilibrium, />., within the slow-binding regime, were analyzed by fitting 
tu an integrated rate equation. The appropriate integrated rate equation for an A + B . * C system going 
to equilibrium does not appear to have been previously derived. The analysis of the slow-binding curves 
v ictded a K t value of 16 pM. in agreement with that of 16 pM determined in the equilibrium titrations, and 
a bimolecular rate constant of association, A on , of 5 0 x 10 ft Ms 1 The experimentally determined k OTI 
and K t indicate that the dissociation rate constant. k M t, is 78 x 1 0 * s ; (/, , = 1 50 mm). The slow-binding 
curves are shown here to fit a simple H + I . * 1:1 model, indicating that it is not necessary to invoke a two-step 
mechanism to explain the inhibition kinetics 
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FiGURfc 1: Structure of franj- Pro- boroPro showing chiral centers. 
The coupling of L-Pro with race mic i n-boroPro is expected to yield 
a mixture uf l*o diaslereomeri L-Pro-L-boroPro and L-Pro-o 
boroPro. 

(Schonera/., 1985; Flentke *r a/., 1991). Specific inhibitors 
of DP IV arc therefore of some interest, both as tools to help 
elucidate the biological role or roles of DP IV and as potential 
therapeutic agents. 

We have previously reported the synthesis and a preliminary 
kinetic characterization of two potent inhibitors of DP IV, 
Ala-boroPro and Pro-boroPro (boroPro refers to an analog of 
pruiine in which I he car boxy late group is replaced by a boron y! 
group) (Bachovchin ef a/ , 1990; Flentke f/j/., 1991). Thc.ic 
inhibitors have an immunosuppressant activity, suppressing 
antigen-induced T-cell proliferation in T-ccIl culture systems 
( Flentke et al„ 1 99 1 ) and antibody production in mice ( Kubota 
et a/., 1992). These findings lend support to the hypothesis 
that DP IV plays a role in T-cell proliferation and suggest 
that DP IV inhibitors may be of therapeutic value. 

Ala-boroPro and Pro-boroPro belong to a class of serine 
protease inhibitors known as peptide boronic acids (Keltner 
& Shenvi, 1984). Inhibitors of this class can have remarkably 
high affinities for their target enzymes. For example, MeO- 
Suc-Ala-Ala-Pro-boroPhe inhibits chymotrypsin with a K x of 
1 60 pM ( Kettner & Shenvi, 1 984),and Ac-D-Phe-Pro-boroArg 
inhibits thrombin with a K\ of 41 pM (Kettner et a/., 1990). 
The potency of these inhibitors h widely attributed to the 
abi.ity of the boronyl group to form a tetrahedral adduct with 
the active site serine, which closely mimics the transition state 
of the enzyme-catalyzed reaction (Koehler & Lienhard, 1971; 
Ltndquist & Terry, 1974; Rawn & Lienhard, 1974; Philipp 
& Maripuri, 1981; Bachovchin et a/., 1988). The peptide 
moiety, however, must also contribute importantly to the 
affinity, as simple alkyl- and arylboronic acids are many orders 
of magnitude less effective as inhibitors. X-ray crystallography 
and NMR spectroscopy have confirmed the presence of a 
boron-serine tetrahedral adduct in several serine protease 
peptide boronic acid inhibitor complexes. However, NMR 
spectroscopy has also demonstrated that in certain cases 
tetrahedral boron histidine adducts are formed { Bachovchin 
et ai., 1988; Tsilikounas et al., 1992). 

The more potent peptide boronic acid inhibitors usually 
inhibit their target enzymes in a time-dependent manner 
( Kettner & Shenvi, 1 984; Shenvi, 1 986; Kettner et a/., 1 988, 
1990), a phenomenon known as slow-binding inhibition 
(reviewed in Morrison and Walsh ( 1 988)]. Both Ala-boroPro 
and Pro-boroPro are slow-binding inhibitors of DP IV 
Morrison and Walsh (1988) have postulated that most, if not 
all. slow-binding inhibitors bind to their target enzymes in 
two steps. f> , the inhibitor first forms a relatively weak 
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understanding the catalytic mechanism of serine proteases 
and for the rational design of inhibitors. 

A major impediment to the study of slow-binding inhibition 
is that the kinetic analysis is not trivial. The high affinity 
these inhibitors typically have for their target enzymes means 
that kinetic experiments must of»en be carried out under 
conditions where I E and, thus, where the approximation 
that If ree = Iioui. which greatly simplifies the kinetic analysis 
of weaker binding inhibitors, is no longer valid. The time 
dependence of the inhibition further complicates matters 
because it may prevent a steady-state rate from being reached 
until substrate depletion becomes significant. Such a system 
is described by a set of differential equations for which an 
integrated rate equation is not available, although expressions 
have been derived for thecase where substrate depletion is not 
significant during the time course of inhibitor binding (Cha, 
1975, 1976). 

The kinetic analysisof Ala-boroProand Pro-boroPro binding 
to DP IV is even more complicated because these inhibitors 
are unstable, having half-) ives of about 5 and 30 min, 
respectively, at neutral pH . In the preliminary kinetic analysis 
we reported K y values of 2 and 3 nM, respectively, for Ala- 
boroPro and Pro-boroPro, realizing that tnese values sub- 
stantially overestimated the true K\ values owing to the 
simplified way in which K t determinations were carried out 
and to the instability and slow binding of these inhibitors 
( Flentke et al., 1991). The original analyses were also carried 
out with inhibitors which were diastereomeric mixtures, i.e., 
L- Ala-DL-boroPro and L-Pro-DL-boroPro. The expectation is 
that only one of the isomers, presumably the l-l isomer, is the 
active inhibitor. 

Because the potency of these inhibitors is unusually high 
for such small molecules, and because DP IV appears to have 
important biological functions, a more detailed analysis of 
how these small dipeptide boronic acids interact with DP IV 
should be of considerable interest. Here we report (i) the 
purification of l- Pro- L- boroPro ?nd L- Pro- D- boroPro from 
the l-dl diastereomeric mixture an * (ii) a more detailed kinetic 
analysisof each isomer's inhibition of DP IV. To circumvent 
the difficulties outlined above, wr have developed a method 
which exploits the fact that dissociation of the inhibitor from 
the enzyme is a relatively slow process. The method involves 
incubating the enzyme with inhibitor in the absence of 
substrate. The amount of free enzyme at any time can then 
be determined by adding substrate and monitoring the time 
course of the enzyme-catalyzed reaction for a short period 
during which the inhibitor does not have time to measurably 
dissociate from the enzyme These simplified experimental 
conditions allow the derivation of expressions which can be 
used to analyze inhibitor binding under both equilibrium and 
nonequilibrium conditions. Equilibrium conditions here refer 
to experiments in which the preincubation time was sufficient 
for en?yme and inhibitor to reach equilibrium prior to the 
addition of substrate and enzyme assay. Nonequilibrium 
conditions refer to experiments in which the preincubation 
time was insufficient for equilibrium to be reached, and thus 
the system is within the slow-binding time domain. The 
integrated rate equation needed to analyze the nonequilibrium 
data does not appear to have been previously delved and is 
therefore derived here for the first time. This approach and 
the derived equations should prove useful in the analysis of 
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previously (Bachovchin et at., 1990). Analytical and semi- 
preparative CI 8 HPLC were performed on a 250 x 4.6 mm 
5-Mm Nucleosil C18 HPLC column (Alltech Associates Inc., 
Deerfield, IL) using a Hewlett-Packard 1 050quaternary pump 
HPLC equipped with a multiple wavelength detector (Hewlett- 
Packard, Rockville, MD). Several milligrams of the purified 
components could be prepared by repeatedly injecting 0.5 mg 
of the mixture on this column and then pooling and lyophilizing 
the appropriate fractions. The resulting material was redis- 
solvedinO.Ol N HC1. Analytical C 1 8 HPLCchromatograms 
of the purified products are shown in Figure 2. The absolute 
configurations were assigned on the basis of a detailed NMR 
study (J. L. Sudmeicr, W. G. Gutheil, and W. W. Bachovchin, 
unpublished results). An attempt to scale up this purification 
procedure on 200 x 1 0 mm and 400 x 1 0 mm Absorbosphere 
CI 8 HPLC columns (Alttech Associates) did not provide as 
pure a final product. 

Quantitation of Pro-boroPro by Amino Acid Analysis. 
Amino acid analysis was performed by the PITC method 
(Bidlingmeyer et at., 1984). Quantitation was based on 
proline. The boronylprolinc did not appear in this analysis. 

Purification of Pig Kidney DP IV. Pig kidney DP IV was 
prepared as described previously (Wolf et a/., 1978). The 
concentration of DP IV active sites was assessed by stoichi- 
ometric titration with L-Pro-L-boroPro, as described further 
below. 

Standard DP IV Emyme Assays. Standard activity assays 
were performed in 50 mM sodium phosphate (pH 7.5) at 25 
°C with the chromogenic substrate Ala-Pro-p-nitroanilide 
(APPNA)(Bachem Inc., Torrance, CA), monitoring the A A \q 
on a Hewlett-Packard UV-vis spectrometer. The value At - 
8800 M 1 cm' 1 upon hydrolysis of substrate was used to 
calculate rates and concentrations (Erlanger et at., 1961). 
The hydrolysis time course was monitored for 2 min. The 
initial substrate concentration was 73.7 ^M, 5 times the K m 
{vide infra). 

Equilibrium Titrations of DP IV with L-L and L-D Pro- 
boroPro. These experiments were performed by f irst preparing 
a DP IV stock in the assay buffer. The amount of DP IV used 
in each assay was the minimal amount necessary to obtain a 
sufficient absorbance change in 2 min with the substrate for 
accurate quantitation. A series of Pro-boroPro dilutions and 
a blank were prepared in 0.0 1 N HC1. To 0,980 mL of the 
stock-diluted DP IV was added 10 mL of diluted Pro-boroPro, 
the mixture was incubated for 30 min at 25 °C, and the free 
enzyme was assayed by the addition of APPNA in 10 of 
DMF 

Kinetics o f DP IV and L- Pro-boroPro Association. For 
the association kinetics a fluorometric assay with Ala-Pro- 
7-amino-4-(trifiuoromethyl)coumarin (APAFC) (Enzyme 
Systems Products, Livermore, CA) was used. Fluorescence 
was monitored on a Perkin-EImer LS-5 fluorescence spec- 
trometer (Oak Brook, IL) with an excitation wavelength of 
400 nm and a detection wavelength of 505 nm. The response 
was calibrated with 1 ^M 7-amino-4-(trifluoromethyl)- 
coumarin. The experiments were performed by incubating 
DP IV with inhibitor in the absence of substrate. After an 
appropriate time interval, APAFC was added to assay for 
free DP IV. The slow apparent rate of DP IV-inhibitor 
association jr.d"r these dilute conditions gave a time course 

■ '-v - .,!!. - - t m ,,-. > \ kMrr ' , Spci f ! . \rr 

used in the equilibrium binding experiment \ oi.nu i.svi , 
on the diluted DH IV was pcrlornied bv mixing u vvu mi. m 



the diluted enzyme with 10 nL of 1 mM APAFC in DMF in 
a cuvette (10 mM final concentration) and monitoring the 
fluorescence change for 2 min at 25 °C. This value was 
considered to be / = 0 for the binding time course. An aliquot 
of L-Pro-L-boroPro was then added to the remaining 9.010 
mL of diluted DP IV, and 0.990-mL aliquots were removed 
and assayed at intervals with APAFC as above. 

Stability o/l-l and l-d Pro-boroPro. In 0.01 N HCI these 
compounds appeared stable for at least 1 month at room 
temperature. Both the l-l and l-d Pro-boroPro preparations 
lose their DP IV inhibitory activity in the assay buffer at pH 
7.5. To partially characterize this behavior, the time course 
of inactivation was monitored at three different concentra- 
tions; two at relatively low inhibitor concentrations, using 
DP IV inhibition as an indicator of residual inhibitor 
concentration, and one at relatively high inhibitor concen- 
tration, using CI 8 HPLC to determine residual inhibitor 
concentration. For inactivation in the range of inhibitor used 
in the assays above, the inhibitor was diluted into assay buffer 
at a concentration sufficient to give roughly 90% inhibition 
initially (approximately 1 nM for the l-l and 100 nM for the 
l-d), At various time intervals DP IV was added, and after 
a i 5-min incubation, APPNA was added to assay for the free 
enzyme analogous to the procedure used in the equilibrium 
titrations. In a second experiment the inhibitor concentration 
in the pH 7,5 buffer was 100 times the assay concentration. 
At time intervals between 1 and 150 min, 10 nL of the Pro- 
boroPro solution was added to 0.980 uL of assay buffer 
containing DP IV, and this mixture was allowed to incubate 
for 15 min. Substrate was then added to assay for free enzyme. 
The inhibition observed in these experiments was converted 
to the amount of active inhibitor using the inverse to the 
equilibrium relationships as derived in the Theory smtion. 
The inactivation of these compounds was also observed directly 
at higher concentrations (250 uM> by analytical C 1 8 HPLC. 
Half-lives {t\ /2 ) for degradation were determined empirically 
from plotted degradation time courses as the time at which 
one-half of the inhibitor remained. 

Progress Curves for Enzyme + Substrate + Inhibitor 
Assays. These experiments were performed in two ways. One 
was for the enzyme to be added to an assay mixture containing 
a known amount of both the substrate and the inhibitor. The 
other was to incubate the enzyme with inhibitor for 1 5 min 
to establish equilibrium and then add substrate. The control 
for this experiment was to monitor a complete time course for 
the hydrolysis of substrate by enzyme. 

Numerical Inetegration of Rate Equations. Numerical 
integrations of rate equations were performed with the G EAR 
software package (Stabler &. Chesick, 1978; McKinney & 
Weigert, 1986). 

Data Analysis. Data were analyzed by fitting to the 
appropriate equation by derivative-free nonlinear regression 
using the IBM PC based version of the BMDP program AR 
(BMDP Statistical Software. Los Angeles, CA). The equa- 
tions used are derived in the Theory section. 

THEORY 

Derivation of Equations Describing Simple A + B ^ C 
Equilibrium. In the enzymalically monitored equilibrium 
titrations, we are titrating DP IV with Pro-boroPro. The 
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exactly using the quadratic equation. Using Ej to represent 
the total DP IV concentration and /? to represent the total 
inhibitor concentration, the following equations can be derived: 



(£fi 



(£ T + / T + K) - \/((f T + / T + /0 2 -4£ T / T ) 



£ = £ T - (£/) 



(1) 
(2) 
(3) 



The observable is the rate of APPNA hydrolysis, which is 
proportional to £: 



rate = £(SA) 



(4) 



where SA is the specific activity in units of AOD 4 io/min/pM 
DP IV active sites at 73.7 uM APPNA. The experimentally 
variable parameter is h. The adjustable parameters to be fit 
arc Ku Ej (in terms of active sites), and SA. In the case of 
L-PrO'D-boroPro, the less potent inhibitor, the observed binding 
appeared to be due to contamination with the l-l diastereomer. 
This siiuuiion was analyzed by including another adjustable 
parameter, in these equations. The actual amount of 
L-l present therefore was 



/ t U-l) = / t (l-d)%l-l/100 



(5) 



where /t(l-l) is the true total l-l concentration, /t(l-d) is 
the total inhibitor concentration (based upon amino acid 
analysis), and %L-L is an adjustable parameter describing the 
% contamination of l-l in the l-d preparation. 

Inversion of the Equilibrium Equation To Measure the 
Rate of Inactivation of l-Pro-l-boroPro. In the preceding 
section, an equation was derived describing the observed rate 
of DP IV catalyzed substrate turnover as a function of the 
independent variable Ml-l) and the parameters /f„ £ T , and 
SA, which are to be fit to experimental data. Once these 
parameters have been fit, it is possible to find the inverse 
relationship to this equation and, with the fitted parameters, 
to calculate /t(l-l) from an experimentally measured rate as 
required for analysis of the DP IV monitored inactivation of 
Pro-boroPro described above. The following equation is easily 
derived: 

/ t (l-l) = E T - AT, - rate/SA + £ T tf,(SA)/rate (6) 
where the parameters are as defined above 

Derivation of an Integrated Rate Equation for the A + B 
— ■ C System. Surprisingly, the appropriate form of the 
integrated rate equation for this system was not found in a 
number of standard sources. The system 



E + / ^ {El) 



(7) 



is described by the following set of differential equations 
d£/d/ = -* on £/+* ofr (£Y) (8a) 
d//df = iJ7 4 k 0 „(EI) (8b) 
dfFn/dr = k r FI k n {FP> (8c) 
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into the expression for d(£/)/df and then expansion and 
rearrangement gives 

d(£0/d/ = k w {EI) 2 - (* on £ T + * on / T + k m Ki){El) + 

*oo£t't d2) 

This can be rearranged for integration as 
d(£/) 



dt = 



U^ 2 - {k^Ej + kjj + k^iEf) + k^Ejlj 

(13) 

The left side of this differential equation is trivial to integrate. 
The right side is given in standard math tables [CRC 
Handbook of Chemistry and Physics, Vol. 67, p A- 26, eq 
1 10, second equation). (Note that we use q for -q and that 
q in our nomenclature can be shown to always be greater than 
or equal to 0, a prerequisite for using this equation.) 

J Vq L 2cx + b + Vq J 

where x = (EI), X = a + bx + cx 2 , a = k^Ejh, b = -(k on E T 
+ k on l-i + k on Ki), c = *o„, and q = b 2 - 4ac. The appropriate 
integrated expression is therefore (with an initial condition at 
r = 0 of x = 0 (i.e., (£/) = 0)) 

Vq L(2cjc + *>+ Vq)(b- Vq) J 
Rearrangement to solve for x in terms of t gives 

(1 -e l{q)i,1 )(b + Vg) 



2c(e'^''-((6 + v^)/(A-v^))) 



(16) 



To check this result, note that at / = 0, x = 0 as expected. Also 
note that 



~* 2c(») 
-b-Vq 



2c 



(17) 



(18) 



which is equivalent to the equilibrium expression derived above 
(eq 1), also as expected. 

Analysis of Enzyme + Substrate Progress Curves. The 
data were collected at 5-s intervals over the time course of 
these experiments, up to 1 .5 h. Several approaches have been 
described for the analysis of data of this type. Direct fitting 
to the integrated Michaciis Menten equation (Kellershohn 
& Larent, 1985; Cox & Boeker, 1987) is complicated by the 
fact that this equation is a mixture of linear and transcendental 
functions in the dependent variable, and therefore f (time) 
must be fit as a function of P (product concentration). A 
more direct approach is to determine the rate (dP/dt) from 
the data and to fit this directly to the Michaelis-Menten 
equation (Canela & Franco, 1986). We use this approach 
here, but have not found it necessary to use a complicated 
weighting scheme nor to fit the time course data to a polynomial 

>j;Hinn v rxfM< ' derivatives Instead the data in terms of 
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differences using the formula 

(d/'/d/), = - /> M )/0 |+1 -/ M ) (19) 

Parameters (K m and /c« t ) were then determined by fitting to 
the Michaelis-Mcnten equation 

Wm^k^Ej/^ + KJ (20) 

The predicted time course was calculated from the fit 
parameters by summation using the following equations ' At 
= 5s): 
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S Q = S T 



(21) 
(22) 

(23) 

Si'So-P, (24) 
which can also be performed easily in a Lotus 1 23 spreadsheet. 

Altcratively, the predicted time course can be obtained by 
numerical integration. The system 



E + S ^ (ES) -+E + P 



(25) 



is described by the following set of differential equations: 
d£/df = -£S*, + (£5)(*_, + A w ) (26a) 
dS/dt = -ESk { + (ES)k^ (26b) 
d(£S)/dr = £S*, - (£S)<*_, + k m ) (26c) 
dP/dt = (£5)^ t (26d) 

Given the initial values for the concentrations of the com- 
ponents in this system and values for the rate constants, the 
GEAR program will provide a simulated time course for 
comparison with the experimentally determined time course. 

Numerical Simulation of the Model Shown in Figure 7 A. 
The model shown in Figure 7A is described by the following 
set of differential equations: 

d£/d, = -K*, + <£$)<*_, + kM) _ KnEI + koff{E/) 

(27a) 

dS/d/= -ES*, + (£S)* , (27b) 

d(£5)/d/ = ESk t - {ES)[k , + k w ) (27c) 

dP/dt = (ES)k al (27d) 

d//d/= k on EI + k off (Er) (21c) 

d(£/)/d/ = ^£7 * orr (£/) (270 

Note that eqs 27 = eqs 8 + eqs 26, i.e., model(Figure 7A) = 
model(£ + / ^ EI) + model(£ + S ES — £ + P). Also 
note that for the one species common to both models, £, the 
differential eq 27a is obtained as 

d£ / d 'modd(Fi l .r«7A) = d£ V d + /^£/) + 

d£ / d 'mod c1 <f + 5^;5 (28) 

Given values for all of the initial concentration and the rate 
constants describing Figure 7A. it is possible to simulate 

I- si 1 p 

of Pro-boroPro diastercomcrs by C 1 8 H Pi C is shown in Figure 
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Figure 2: C18 HPLC chromatograma showing resolution of Pro- 
boroPro diastereomers. The upper chromatogram is of the starting 
mixture, the middle chromatogram is of the purified L-Ddiastereomer 
and the lower chromatogram is of the purified l-l diastereomer' 
H PLC conditions: solvent A 0.1% trifiuoroacetic acid (TFA) i n H 2 0- 
solvent B 70% acetonitrile/30% HjO/0.086% TFA; gradient 0-2 
mm 0% B, 2-32 min 0-100% B. Only the first 15 min of each 
chromatogram are shown. The base-line disturbance at 8.1 5 min is 
the gradient entering the detector. 
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Figure 3: DP IV monitored inactivation kinetics of l-l (■) and l-d 
(•) Pro-boroPro at the very low concentrations ( ~ 1 nM of l-l) and 
pH (7.5) used in the equilibrium and kinetic experiments. The l-l 
concentration was calculated from the raw data using the inverse of 
the equilibrium equation derived in the Theory section (eq 6). 

2. Attempts to scale up this separation were unsuccessful 
owing to the decreased resolution with larger column diam- 
eters, which could not be overcome by extending the length 
of the column. For separation on analytical columns, about 
0.5 mg of the mixture was loaded for each run. A purity of 
>98% for the purified products was indicated by analytical 
HPLC (Figure 2V 

Stability of Pro-boroPro Diastereomers. The concentration 
dependence of the inactivation of Pro-boroPro was examined 
at very low (~ 1 nM) and low (-100 nM) concentrations of 
Pro-boroPro, at pH 7.5, using DP IV inhibition to monitor the 
residual active inhibitor concentration, and at a relatively high 
( - 250mM) concentration of the inhibitors using C 1 8 HPLC. 
The time courses of the degradation kinetics had the same 
shapes as the curves shown in Figure 3 in all cases, but with 
somewhat different half-lives The inactivation reaction 

■' t '" f -MM,.:; Jin Mis; H ik i * I HC U. 
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;L-0] <nM) 

Figure 4: Equilibrium titration data and fits obtained for l-l and 
L-D Pro-boroPro inhibition of DP IV. Each point for the l-l titration 
is the average of three determinations. Each point for the l-d titration 
is the average of five determinations. (A) Fit of the l-l titration data 
to a simple equilibrium model (eq 4). (B) Fit of the l-d titration data 
to a simple equilibrium model (eq 4). (C) Fit of the l-d titration 
data to a simple equilibrium model (eq 4) where observed inhibition 
by l-d is due to contamination by a fractional amount, %l-l, of l-l 
(eq 5). In this figure, • is used to represent data points. ■ is used 
to represent the residual value between the calculated value based 
upon the fitted parameter and the experimentally observed values, 
and the lines represent the calculated value of the observable based 
upon the fitted parameter values 

measured by HPLC. The inactive material can be reactivaied 
by acidification, a process which is relatively slow (85% in 18 
h at pH 2.0). The inactive material the cyclic structure in 
which the N-terminal nitrogen atom foi ns a covalent bond 
with the boron atom of the boron yl group (J. L. Sudmeier, 
W. G. Gutheil, and W. W. Bachovchin, unpublished results). 

Equilibrium Binding of L-L and L-D Pro-boroPro to DP IV. 
Under the conditions used for the equilibrium binding assays, 
>95% of the final equilibrium inhibition was observed in 3 
min at the lowest inhibitor concentration used. Free enzyme 
was assayed after incubating enzyme and inhibitor for 30 

upward curvature *a> ooscrveu *nicn -aouiu iridic. 

■ iDltnr *, is OCinf! Oispiaceu *H MjftSir,!:; r; \r.:s -trm- n^.:: 

'■ he data obtained trom trie equilibrium Dinding assays anu 
the statistically fit parameters are shown in Figure 4 I hese 
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Figure 5: Kinetics of L-Pro-L-boroPro binding to DP IV at two 
inhibitor concentrations. Inhibitor concentrations for the two 
experiments are shown in the inset. Date were fit toeq 17. Calculated 
*orr ■= - 78 X lfr* r". 

experiments yield a Ki for L-Pro-L-boroPro of 16 pM (Figure 
4A). I his value is much icss ihan the estimated DP IV 
concentration of 570 pM (SE « 20) (in terms of binding sites) 
and establishes that the titration was well into the tight binding 
regime. The L-Pro-D-boroPro binding data fit poorly to the 
simple equilibrium model ( Figure 4B) but very well to a model 
in which the observed inhibition is assigned to small amounts 
of contaminating L-Pro-L-boroPro( Figure 4C). Thisanalysis 
indicated that our preparation of L-Pro-r>boroPro contained 
0.59% of the l-l diastereomer and that the inhibition observed 
with this preparation is due to the contaminating l-l 
diastereomer. This level of contamination was confirmed by 
analytical HPLC. 

Association and Dissociation Rate Constants for L-Pro- 
L-boroPro Binding to DP I V. To determine values for k M and 
£ 0 ff, a more sensitive fluorescent assay was employed to allow 
the time course of inhibitor binding to DP IV to be monitored 
under conditions of low DP IV concentration ( 1 1 .4 pM active 
sites). A fit of these data to the integrated rate equation 
derived above (Figure 5) yielded a bimolecular association 
rate constant, of 5.0 X 10* M" 1 s" 1 (SE = 0.2), a K> of 
16 pM (SE = 7) (a value identical with that determined by 
equilibrium titration), and a unimoIecuLr off rate, k^ 
calculated from the *„, and K» of 78 X lO^s 1 . Thecaiculated 
off rate indicates that the t i/2 for dissociation is 150 min, 
which is slow relative to the 2 min required for the enzymatic 
assays employed in these experiments, thereby confirming 
the original assumption on which this experimental approach 
was based. This t\ ft for dissociation is much longer than the 
t ] /2 for the inactivation of free inhibitor and therefore indicates 
that inhibitor in the enzyme-inhibitor c iplex is more stable 
than the free inhibitor. 

Kinetics of APPNA Hydrolysis by DP I V. Progress curves 
for the hydrolysis of APPNA starting at the standard assay 
concentration of 73.7 APPNA are shown in Figure 6A. 
The fit to the Michaelis-Menten equation using the method 
described above is very good (Figure 6B), and only a slight 
difference between the experimental and fit curves is dis- 
cernible. Product inhibition therefore appears negligible in 
this case This analysis yielded a k ax of 90 8 s 1 (SE = 0.9) 

r i^rt^ , un f '■ ■ i. n:\ mt upstru, urnm 

,n\tina \ umem ai Simulation > . » simulate inr fivOrolvsi 
; substrate and also the binding and dissociation ot inhibitor 
in trie presence of substrate requires the extraction ot the 
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Figure 6: Progress curve for the hydrolysis of A PPN A in the absence 
and presence of L-Pro-L-boroPro and fitted and simulated progress 
curves. (A) Progress curve data for hydrolysis of APPNA by DP 
IV. The initial concentration of APPNA was 73.7 mM and that of 
Uf i V was 5"0 pM. (3) Curve resulting from the fit of the progress 
curve to the Michaelis-Menten equation (eq 20) as described in the 
text. The fit parameters were K m m 14.3 mM (SE = 0.4) and *«, 
- 90.8 s- 1 (SE - 0.9). (C) Substrate hydrolysis progress curve 
simulated by numerical integration of eqs 26a-d using k x = 6.34 x 
10* M" 1 s- 1 , * 0, and k m ■ 90.8 s~\ as discussed in the text. The 
curves A-C are essentially indistinguishable. (D) Experimental data 
from the experiment where DP IV is added to a fresh mixture of 
APPNA(73.7 M M)andL-Pro-L-boroPro(4.76nM). (E)Simulation 
of experiment D by numerical integration of eqs 27a-f with * (t k i, 
and * M i as in C and k w = 5.0 X 10* M-' s-' and k^ » 78 x l(H s '. 
(F) Experimental data for the experiment where DP IV is preincubated 
with L-Pro-L-boroPro (4.76 nM) for 1 5 min and then APPNA (73.7 
mM) is added. (G) Simulation of experiment F by numerical 
integration. 

component rate constants from the K m : 



Equation 29 can be rearranged to give 



(29) 



(30) 



Since k-i must be greater than or equal to 0, k } must be 
greater than or equal to 6.34 X 1 0 6 M _I s"' on the basis of the 
values of K m and k<* t obtained above for DP IV and APPNA. 
The minimal model for simulating the kinetics of substrate 
hydrolysis employs values of k.\ - 0, k x = 6.34 X I0 6 M" 1 
s and kctt = 90.8 s 1 . A simulation of the time course nf 
APPNA hydrolysis by DP IV using these values by numerical 
integration of eqs 26 gives the result shown in Figure 6, curve 
C. It is clear that the parameter values obtained using these 
methods accurately account for the DP IV catalyzed progress 
curves for APPNA hydrolysis. These values for k. u and 
*<ui, together with the determined k M and k otf values for l-l 
Pro-boroPro binding to DP IV, allow the kinetics of a ternary 
enzyme + substrate + inhibitor experiment, or of substrate 
being added to preequi libra ted enzyme + inhibitor, to be 
simulated by numerical integration of the differential equations 
describing Figure 7A (eqs 27). The results from these 
experiments and simulations are shown in Figure 6, curves 

uratm r etiCL: me simplified experiments from *mch tru- 
c 'tni.ii net]. ;r:c, u. : >,,n a^Lratei\ predict trie more complt 
'i/vme + substrate + inhibitor experiments It should De- 
noted that in our simulations the depletion of free enzyme bv 
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Figure 7: Two models for inhibition which have been proposed and 
used for the analysis of slow-binding inhibitor data: (A) simple model; 
(B) complex model. 

the formation of the ES complex is taken into account and 
that formation of the ES complex is rapid and reaches a steady- 
state concentration in about 10 ms. Increasing the rate at 
which ES reaches the steady sta f * by increasing the value of 
k. u and the value of k x to satisfy eq 30, has no significant 
effect on either the simulated substrate hydrolysis progress 
curve or the simulated enzyme + substrate + inhibitor progress 
curves. 

DISCUSSION 

The slow rate at which L-Pro-L-boroPro dissociates from 
DP IV allows the amount of free enzyme to be sampled during 
the course of an enzyme plus inhibitor preincubation by the 
addition of substrate, without the added substrate perturbing 
the position of the enzyme-inhibitor equilibrium, over the 
coui e nf a 2-min assay. This simplification in the experimental 
desigr for a relatively simple mathematical treatment 

of inhibitor binding to DP IV as an A + B ^ C system at 
equilibrium or approaching equilibrium. Analytical equations 
for both situations are derived. The integrated rate equation 
describing the kinetics of :uch a simple system approaching 
equilibrium did not appear in standard sources, although we 
did not search extensively into the older literature. The 
experimental data were fit directly to these equations, and the 
reliability of the values for the fit parameters was assessed 
from the standard errors and from the correlation matrix from 
the analysis. Other parameters, such as stoichiometry or 
contamination level, can be included in such an analysis 
relatively easily, as demonstrated above for the determination 
of the active site concentration and for the determination that 
the observed inhibition by our preparation of L-Pro-D-boroPro 
is due to contaminating L-Pro-L-boroPro. 

The experiments were designed to minimize the effect of 
inhibitor instability on the acquired data, since it was known 
from previous studies that Pro-boroPro had a half-life of 
between 30 min and I hatpH values of7.0 and above (Flentke 
rt at, 1991) The equilbrium titrations were therefore 

ninitor iu ensure thai association is complete in 1 nun i n ^ 
"UUiithr nun experiments should have suifereti nine '■>? 
perturbation from inhibitor instability On the other hand, 
the kinetic experiments had to be conducted at enzyme and 
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inhibitor concentrations low enough to allow the rate of 
association to be measured and thus were more susceptible to 
the effect of inhibitor instability. The effect of inhibitor 
instability was minimized by restricting the time course to 30 
min. during which -20% of the inhibitor would have 
decomposed (Figure 3). Thus, even in the kinetic experiments 
the effects from inhibitor instability should have been small. 

The results from the application of this approach to the 
binding of Pro-boroPro to DP IV demonstrate that L-Pro-L- 
boroPro is a tight binding inhibitor of DP IV with a K\ of 16 
pM . This value is substantially lower than the original estimate 
of 3 nM, as expected, and more accurately reflects the true 
affinity of L-Pro-L-boroPro for DP IV. The results also 
demonstrate that the l-d diastcreomer is a much weaker 
inhibitor than the i-L, if it is an inhibitor at all, as the analysis 
of the data demonstrates that all of the observed inhibition 
can be attributed to the presence of ^ 0.59% of the l-l 
diastereomer ( Figure 4B,C). The bimolecular for L-Pro- 
L-boroPro of 5.0 x |0 6 M 1 s 1 , obtained from the kinetic 
experiments, is close to the values of (5-9) X 10 6 M' 1 s _l 
measured for peptidyl boroArg based inhibitors of thrombin 
(Kcttr.cr c: a!., 1990). The k^ cf 78 X !0~* s" 1 , calculated 
from the experimentally determined k 0R and K u indicates that 
the tight binding of these inhibitors is a consequence of their 
slow dissociation rates. The slow dissociation of inhibitor (tij 2 
= 150 min) validates the experimental approach used in this 
study of preincubating enzyme and inhibitor followed by 
assaying the enzyme by the addition of substrate. 

The tight binding behavior of L-Pro-L-boroPro with DP IV 
allows it to be used as an active site titrant. The accuracy 
with which we can measure the concentration of DP IV active 
sites should be limited primarily by the accuracy with which 
wecan measure theconcentration of thestock L-Pro-L-boroPro 
solution. The concentration of this stock solution was 
determined from three replicate analyses to be 4.76 mM with 
a standard error of 0.05 or 1.1%. The results of an active site 
titration of DP IV with L-Pro-L-boroPro together with the 
published A< of 8800 M 1 cm 1 at 410 nm for substrate 
hydrolysis (Erlanger et a/., 1961), and the progress curve 
analyses shown in Figure 6B, yield a or turnover number 
of 90.8 s" 1 . This is significantly higher than the 54.6 s 1 
previously reported for this enzyme-substrate pair obtained 
underslightly different experimental conditions (Heins etal. y 
1988) The values obtained in this study are less subject to 
errors owing to protein impurities or denaturation than values 
based upon protein concentration measurements. Because 
[ - Pro-L-boroPro is a transition-state (or intermediate) analog, 
only catalyiica!!) active enzyme is expected to hind the 
inhibitor with high affinity. We cannot, however, rule out 
the possibility that species other than calalytically active DP 
IV might bind to the inhibitor with sufficient affinity to 
introduce some error into the use of the stoichiometry for 
inhibitor binding as a measure of calalytically active DP IV 
The estimate of the value for * Mt of 90.8 s 1 for APPNA 
hydrolysis under these conditions should therefore be con- 
sidered as the best current estimate. The true value will in 
fret be higher if species other than catr^'ytically active DP IV 
bind this inhibitor with high affinity The K m of 14.3 jiM 
obtained here is similar to the previously reported value of 
16 6 ,;M (Hems et al , 1988) 
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formed which then undergoes slow conversion to a tighter 
complex [Figure 7B, mechanism B of Morrison and Walsh 
( 1 988)]. The decision of which mechanism is followed should, 
in principle, be a simple matter of determining which 
mechanism better accounts for the experimentally observed 
progress curves. In practice, the discrimination between the 
one- and two-step mechanisms is not so simple. A major 
problem is that the integrated forms of the differential 
equations describing the two-step mechanism (mechanism B) 
have not been derived owing to their complexity. Conse- 
quently, alternative methods have been sought, and the most 
often used method has been to fit the observed slow-binding 
progress curves (of enzyme + substrate + inhibitor) to the 
integrated form of the one-step model (Cha. 1976; Williams 
et al., 1979). This approach yields values for the initial rate 
at the beginning of the progress curve, »u the final steady- 
state rate, v %% and the apparent first-order rate constant for 
the approach to steady state, k obi , A dependence of e 0 
(Williams et al,. 1979) or a nonlinear dependence of * obi on 
inhibitor concentration is taken as evidence for a two-step 
mechanism (mechanism B). Williams et al. (1979) have 
carried out a more rigorous analysis of the slow-binding 
inhibi# ion of dihydrofolate reductase by methotrexate by fitting 
the data to a simulation generated by numerical integration. 
Such an approach has thus far been reported only once. 
Williams et a ' (1 979) concluded that the observed inhibition 
kinetics of DH FR (dihydrofolate reductase) with methotrexate 
was consistent with a two-step mechanism, but they did not 
show that the fit to a two-step mechanism was statistically 
better than that to a one-step mechanism. 

The approach we have employed here of preincubating DP 
IV and L-Pro-L-boroPro in the absence of substrate yields 
slow-binding curves that both qualitatively and quantitatively 
resemble slow-binding curves generated in the presence of 
substrate. A major advantage of this approach is that, owing 
to the absence of substrate, mathematical analysis of the slow- 
binding inhibition becomes greatly simplified. The integrated 
equation for a one-step mechanism in the absence of substrate, 
however, does not appear to have been derived before and is 
given in eq 1 6. A fit of the observed slow-binding kinetics in 
the absence of substrate to eq 16 is quite good (Figure 5) and 
yields a K) value which agrees with that obtained in the 
equilibrium incubation experiments. This analysis thus 
indicates that l-l Pro-boroPro binds to DP IV via a one-step 
mechanism and that the time dependence of the inhibition 
arises simply as a consequence of low enzyme and inhibitor 
concentrations. This conclusion is at odds with previous studies 
of slow-binding inhibition of serine proteases by peptide boronic 
acids, which have generally favored a two-step mechanism 
(Kcttncr & Shenvi, !984;Shenvi. 1 986; Kettner et al , 19^8, 
1990). Intuitively, it seei . that the slow-binding inhibition 
observed for this class of inhibitors should have a common 
mechanism and, therefore, that perhaps some of the earlier 
kinetics supporting two-step binding should be reexamined 
using the methods reported here. Such studies are in progress 

The parameters obtained for the analysis of the DP IV + 
L-Pro-L-boroPro data in the absence of substrate accurately 
account for the data obtained in these experiments (Figures 
4 and 5) A problem arises, however, when these parameters 
:ire used to simulate the inhibition kinetics observed for the 

<' - h ■ ,t m \a i i v ■. . :t r-' 'Mi.t i : : c \ ;><:: : 1 1 l n i - 1 ii \ >i ti.se r \. :. 
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holding the substrate parameters fixed, yields a K x for the 
inhibitor of 510 pM, a value substantially higher than that of 
16 pM obtained with the simplified system. As the inhibitor 
concentration is increased, the apparent K\ value decreases, 
suggesting that the K\ value obtained with the simplified 
experimental approach represents a limiting value which would 
eventually be realized at a sufficiently high inhibitor con- 
centration. The apparent dependence of K x on the inhibitor 
concentration is inconsistent with either a one- or two-step 
inhibitor binding mechanism and must therefore reflect some 
as yet unidentified complication in the enzyme + substrate 
+ inhibitor mixtures which is obviated in the preincubation 
experiments. The possibility that a component of the substrate 
solution accelerated the inactivation of the inhibitor has been 
ruled out (data not shown). Other possibilities include some 
typeof allosteric interaction between the two identical subunits 
of DP IV or hysteresis. 

The experimental approach reported here can be applied 
to enzyme-inhibitor systems in which the dissociation of the 
inhibitor is slow relative to the time required to assay the free 
enzyme. Such an approach has advantages and disadvantages 
over the traditional method. The main disadvantage is that 
this approach requires the / t/2 for dissociation of the inhibitor 
to be slow relative to the time necessary for measurement of 
the enzyme activity. A second disadvantage is that the ki netic 
measurements require a very sensitive assay capable of 
measuring the small amounts of enzyme present under kinetic 
conditions. The major advantage of this approach is that it 
allows data obtained under both equilibrium and kinetic 
conditions to be fit directly to analytically derived equations 
formulated in terms of fundamental physical constants. 
Another advantage is that it yields kinetic parameters for 
inhibitor binding under conditions where complications arising 
from substrate are eliminated. Such an analysis provides a 
useful starting point for dissecting more complex kinetic 
systems such as the one described here. We expect that, in 
the absence of complicating factors, this approach will provide 
parameters which fully account for the behavior of the enzyme 
+ substrate + inhibitor system. The numerical simulation of 
complex kinetic systems using experimentally derived rate 
and equilibrium constants is not commonly used to verify that 
such constants accurately reflect the data from which they 
were obtained. The results reported here demonstrate that 
such an exercise may reveal more complicated interactions 
between components than originally anticipated. 
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